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A search for the production of neutral Higgs bosons $ decaying into t + t~ final states in pp 
collisions at a center-of-mass energy of 1.96 TeV is presented. The data, corresponding to an 
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integrated luminosity of approximately 325 pb" 1 , were collected by the D0 experiment at the 
Fermilab Tevatron Collider. Since no excess compared to the expectation from standard model 
processes is found, limits on the production cross section times branching ratio are set. The results 
are combined with those obtained from the D0 search for §b(b) — > 666(6) and are interpreted in the 
minimal supersymmetric standard model. 

PACS numbers: 14.80.Bn, 14.80.Cp, 14.60.Fg, 13.85.Rm, 12.60.Fr, 12.60.Jv 



Final states leading to high-mass tau lepton pairs can 
arise from various physics processes beyond the standard 
model (SM) including the production of neutral Higgs 
bosons (generally denoted as <P). Higgs bosons are an 
essential ingredient of electroweak symmetry breaking in 
the SM, but so far remain unobserved experimentally. 
A search for Higgs bosons decaying to tau leptons is of 
particular interest in models with more than one Higgs 
doublet, where production rates for pp — > <P — ► rr can 
potentially be large enough for an observation at the Fer- 
milab Tevatron Collider. For instance, the minimal su- 
persymmetric standard model (MSSM) [l| contains two 
complex Higgs doublets, leading to two neutral CP-even 
(h,H), one CP-odd (A), and a pair of charged (-H^) 
Higgs bosons. At tree level, the Higgs sector of the 
MSSM is fully specified by two parameters, generally cho- 
sen to be Ma, the mass of the CP-odd Higgs boson, and 
tan /?, the ratio of the vacuum expectation values of the 
two Higgs doublets. At large tan/?, the coupling of the 
neutral Higgs bosons to down-type quarks and charged 
leptons is strongly enhanced, leading to sizeable cross 
sections and increased decay rates to the third gener- 
ation tau lepton and bottom quark. MSSM scenarios 
with large tan [3 are of considerable interest since they 
can provide a viable dark matter candidate 0] ■ 

Searches for neutral MSSM Higgs bosons have been 
conducted at LEP 0] and at the Tevatron 0, Q . In this 
Letter a search for $ — > rr decays is presented. At least 
one of the tau leptons is required to decay leptonically, 
leading to final states containing erh, fJ-Th and e/i, where 
Th represents a hadronically decaying tau lepton. 

The data were collected at the Fermilab Tevatron 
Collider between September 2002 and August 2004 at 
y/s = 1.96 TeV and correspond to integrated luminosities 
of 328 pb" 1 , 299 pb" 1 , and 348 pb" 1 for the er h , fir h and 
e/Lt final states, respectively. Final states with two elec- 
trons or two muons have a small signal-to-background 
ratio due to the small branching fraction and the large 
background from Z/j* production, and are therefore not 
considered. 

A thorough description of the D0 detector can be 
found in Ref. Briefly, the detector consists of a 

magnetic central tracking system surrounded by a liquid- 
argon and uranium calorimeter and a toroidal muon spec- 
trometer. The central tracking system comprises a sili- 
con microstrip tracker (SMT) and a central fiber tracker 
(CFT), both located within a 2 T magnetic field provided 
by a superconducting solenoidal magnet. The SMT and 



CFT designs were optimized to provide precise tracking 
and vertexing capabilities over the pseudorapidity range 
\r)\ < 2.5, where r\ = — ln(tan(6'/2)) and 9 is the polar 
angle with respect to the proton beam. The calorime- 
ter is divided into a central section covering < 1.1, 
and two end calorimeters (EC) that extend coverage to 
|?7| w 4.2. A muon system, at | -7-7 1 < 2, consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T toroids, followed by two similar layers af- 
ter the toroids. The luminosity is measured by detect- 
ing inelastic pp scattering processes in plastic scintilla- 
tor arrays located in front of the EC cryostats, covering 
2.7 < M < 4.4. 

The er/, and the fxr^ analyses rely on single electron 
and single muon triggers, respectively, while the e/i anal- 
ysis uses dilepton triggers. Signal and SM processes are 
modeled using the pythia 6.202 Q Monte Carlo (MC) 
generator, followed by a GEANT-based Q simulation of 
the D0 detector geometry. All background processes, 
apart from QCD multijet production, are normalized 
using cross sections calculated at next-to-leading order 
(NLO) and next-to-NLO (for Z boson, W boson, and 
Drell-Yan production) based on the CTEQ5 [9j parton 
distribution functions (PDF). 

The normalization and shape of background contri- 
butions from QCD multijet production, where jets are 
misidentified as leptons, are estimated from the data by 
using like-sign e and candidate events (er^ analysis) 
or by selecting background samples by inverting lepton 
identification criteria (fiTh and efi analyses). These sam- 
ples are normalized to the data at an early stage of the 
selection in a region of phase space dominated by multijet 
production. 

Isolated electrons are reconstructed based on their 
characteristic energy deposition in the calorimeter, in- 
cluding the transverse and longitudinal shower profile. In 
addition, a track must point to the energy deposition in 
the calorimeter, and the track momentum and calorime- 
ter energy must be consistent. Further rejection against 
background from photons and jets is achieved by using 
a likelihood discriminant, which is exploiting character- 
istic calorimeter and tracking information. Muons are 
selected using tracks in the central tracking detector in 
combination with patterns of hits in the muon detector. 
Muons are required to be isolated in both the calorime- 
ter and the tracker. Reconstruction efficiencies for both 
leptons are measured using data. 

A hadronically decaying tau lepton is characterized by 
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a narrow isolated jet with low track multiplicity. The 
tau reconstruction is either seeded by calorimeter energy 
clusters or tracks Three r- types are distinguished: 

• r-type 1: a single track with energy deposition in 
the hadronic calorimeter (1-prong, 7r -like); 

• r-type 2: a single track with energy deposition in 
the hadronic and the electromagnetic calorimeter 
(1-prong, p ± -like); 

• r-type 3: two or three tracks with an invariant mass 
below 1.1 or 1.7 GeV, respectively (3-prong). 

A set of neural networks, one for each r-type, has been 
developed based on further discriminating variables. The 
neural networks were used elsewhere for a cross section 
measurement of the process Z/j* — > rr fioj ]. The input 
variables exploit the differences between hadronically de- 
caying tau leptons and jets in the longitudinal and trans- 
verse shower shape as well as differences in the isolation 
in the calorimeter and the tracker. The training of the 
neural networks is performed using multijet events from 
data as the background sample and tau MC events as 
signal, resulting in a network output close to one for tau 
candidates and close to zero for background. Both train- 
ing samples cover the kinematic region of interest for this 
analysis. For r- types 1 and 2, hadronic tau candidates 
are required to have a neural network output greater than 
0.9. Due to the larger background contamination, this 
cut value is tightened to 0.95 for r-type 3. 

Electrons and muons can be misidentified as one-prong 
hadronic tau decays. Hadronically decaying tau lep- 
tons deposit a significant fraction of their energy in the 
hadronic part of the calorimeter. To reject electrons, 
the ratio between the transverse energy in the hadronic 
calorimeter and the transverse momentum of the tau 
track is required to be larger than 0.4. With a smaller 
rate, background from muons occurs in r-types 1 and 2 
in the /ir^ analysis. This background is suppressed by re- 
jecting tau candidates to which a muon can be matched. 

The signal is characterized by two leptons, missing 
transverse energy and little jet activity. It would 
stand out as an enhancement above the background from 
SM processes in the visible mass 

M vis = ^(P T1 + P T2 + f T ) 2 , (1) 

calculated using the four vectors of the visible tau de- 
cay products P Tl2 an d of the missing momentum fa = 
(^T^XT^yjO). Tfi x and Tf, y indicate the components 
of |?t- For the optimization of the signal selection, 
only the high mass region is used, which is defined as 
Af V i S > 120 GeV in the er^ and fiTh analyses and as 
Af v i s > 110 GeV in the e/i analysis. 

In the eTh and /j,Th analyses, an isolated lepton (e, jj,) 
and an isolated hadronic tau with transverse momenta 
above 14 GeV and 20 GeV respectively are required. In 



TABLE I: Numbers of events observed in data and ex- 
pected for background and the efficiency for a signal with 
M,p = 150 GeV for the three analysis channels, with statisti- 
cal and systematic uncertainties added in quadrature. 



Analysis 


eTh 




e /i 


Data 


337 


575 


41 


QCD 


144 ± 19 


62± 7 


2.1± 0.4 


— » TT 


130 ± 17 


492 ± 53 


39± 5 


Z/7* — > ee, fifi 


12± 2 


5± 1 


0.6± 0.1 


W — ► eu, fiv, tv 


9± 1 


14± 2 


0.3± 0.2 


Di-boson 


0.4± 0.1 


3.1 ± 0.3 


1.0 ± 0.1 


tt 


0.3± 0.1 


1.2 ± 0.2 


0.06 ± 0.02 


Total expected 


296 ± 38 


576 ± 62 


44± 5 


Efficiency % 


3.6± 0.4 


8.6± 0.8 


4.3± 0.5 



addition to the irreducible background from Z/'y* — > tt 
production, a W —* tv decay can be misidentified as a 
high-mass di-tau event if it is produced in association 
with an energetic jet that is misidentified as a hadronic 
tau decay. In these events, a strongly boosted W bo- 
son recoils against the jet, and the mass of the W bo- 
son can be reconstructed in the following approxima- 
tion My/ 1 = \j2 E v E e /^ (1 - cos A(/>), where the az- 
imuthal angle A</> is between the lepton and $t, an d 
E v = E T ■ E e /^/E^. is required to be less than 

20 GeV. 

In the efj, analysis, two isolated leptons each with 
Pt > 14 GeV are required. The dominant background 
contributions after the lepton selection come from the ir- 
reducible Z/j* — > rr process, followed by WW, WZ, tt, 
W — > Iv, and multijet events. In this analysis the multi- 
jet background is suppressed by requiring $t > 14 GeV. 
Background from W+jet events can be reduced using 

the transverse mass M^/ M = \J 2 p^ $t (1 — cos A(j>) 
by requiring that either Mf, < 10 GeV or < 10 GeV. 
Furthermore the minimum angle between the leptons and 
the Et vector, min[A0(e,^r), A0(^,$r)], has to be be 
smaller than 0.3. Finally, contributions from tt back- 
ground are suppressed by a cut on the scalar sum of the 
transverse momenta of all jets in the event Ht < 70 GeV. 

The numbers of events observed in the data and those 
expected from the various SM processes show good agree- 
ment, as can be seen in Table U and Fig. [TJ The estimate 
of the expected numbers of background and signal events 
depends on numerous measurements that introduce a sys- 
tematic uncertainty: integrated luminosity (6.5%), trig- 
ger efficiency (l%-4%), lepton identification and recon- 
struction efficiencies (2%-5%), jet and tau energy calibra- 
tion (2%-6%), PDF uncertainty (3%-4%), and modeling 
of multijet background (2%~9%). All except the last one 
are correlated between the three final states. 

The efficiencies for a Higgs boson signal are found to 
vary between 1.6%, 4.0%, and 1.2% for = 100 GeV 
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FIG. 1: The distribution of the visible mass M v i s for the 
two final states involving hadronic tau decays and for the 
efi final state. The Higgs signal is normalized to the cross 
section excluded by this analysis. The upper distribution 
shows the subsample with the largest signal-to-background 
ratio (M^/ 1 < 6 GeV). The highest bin includes the overflow, 
the indicated luminosity represents the average of the three 
final states. 



and 8.3%, 13.6%, and 9.3% for M<p = 300 GeV for the 
eTh, fiTh, and e/i analyses respectively. Since no signifi- 
cant evidence for the production of neutral Higgs bosons 
with decays <P — ► tt is observed, upper limits on the 
production cross section times branching ratio are ex- 
tracted as a function of M$. In order to maximize the 
sensitivity (expected limit), the event samples of the eTh 
and /zt/j analyses are split into subsamples according to 
different signal-to-background ratios: The subsamples 
are separated by r-type and by Mw (M^ < 6 GeV, 




100 150 200 250 300 

[GeV] 

FIG. 2: The observed and expected 95% CL limits on the 
cross section times branching ratio for & — > tt production 
as a function of M$ assuming a narrow width of the Higgs 
boson. The error bands include systematic and statistical 
uncertainties. CDF curves are taken from [5(, where data 
corresponding to an integrated luminosity of 310 pb^ 1 is used. 



6 < M^/ 1 < 20 GeV). Furthermore the differences in 
shape between signal and background are exploited by 
using the information of the full mass spectrum of M v ; s 
in the limit calculation. Both the expected and the ob- 
served limits on the cross section times branching ratio at 
the 95% confidence level (CL), calculated using the mod- 
ified frequentist approach [11| , are presented in Fig. 

In the MSSM, the masses and couplings of the Higgs 
bosons depend, in addition to tan/3 and Ma, on the 
SUSY parameters through radiative corrections. In a 
constrained model, where unification of the SU(2) and 
U(l) gaugino masses is assumed, the most relevant pa- 
rameters are the mixing parameter X t , the Higgs mass 
parameter \x, the gaugino mass term M2, the gluino 
mass m g , and a common scalar mass Msusy- Lim- 
its on tan/3 as a function of Ma are derived for two 
scenarios assuming a CP-conserving Higgs sector: the 
so-called m™ ax scenario (with the parameters MgusY= 
1 TeV, X t = 2 TeV, M 2 = 0.2 TeV, n = ±0.2 TeV, 
and m g — 0.8 TeV) and the no-mixing scenario (with 
the parameters M SU sy= 2 TeV, X t = 0, M 2 = 0.2 TeV, 
H = ±0.2 TeV, and m g = 1.6 TeV) 0. The production 
cross sections, widths, and branching ratios for the Higgs 
bosons are calculated over the mass range from 90 to 
300 GeV using the feynhiggs program where the 
complete set of one-loop corrections and all known two- 
loop corrections are incorporated. The contributions of 
SUSY particles in the loop of the gluon fusion process are 
taken into account, as well as mass- and tan /3-dependent 
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FIG. 3: Region, which is excluded at 95% CL, in the 
(Ma, tan (3) plane for the m™ ax and the no-mixing scenario 
for ft = +0.2 TeV and fj, = -0.2 TeV (m top = 172.7 GeV). 
The results obtained in the present analysis (labeled as D0 
tt) are combined with those obtained in the <l>b(b) — > bbb(b) 
search The LEP limits [3J have been extrapolated for 
tan/3 > 50. 



decay widths. In the region of large tan/3, the A boson 
is nearly degenerate in mass with either the h or the H 
boson, and their production cross sections are added. 

Fig. [3] shows the D0 results obtained in the present 
analysis in combination with those obtained in the 
<Pb(b) — > bbb(b) search [4j, which are re- interpreted in 
the MSSM scenarios used in this Letter. The combined 



result currently represents the most stringent limit on 
the production of neutral MSSM Higgs bosons at hadron 
colliders. 
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